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Existing lattice data on the QCD phase transition are analyzed in renormalized perturbation 
theory. In quenched QCD it is found that T c scales for lattices with only 3 time slices, and that 
T c /A-jjg- = 1.15±0.05. A preliminary estimate in QCD with two flavours of dynamical quarks shows 
that this ratio depends on the quark mass. For realistic quark masses we estimate T c /Aj^ = 0.49 ± 
0.02. We also investigate the equation of state in quenched QCD at 1-loop order in renormalised 
perturbation theory. 
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The QCD phase transition temperature (T c ) is a fundamental constant of the hadronic world, and will soon be 
accessible to experiments. It has also been the target of many lattice computations. Current practice is to express T c 
in units of the mass of the rho meson (M p ) or the square root of the string tension (^/a) |Q. However, in the last few 
years it has become clear that the renormalised QCD coupling, a s , measured on fairly coarse lattices J2| yield results 
comparable to those obtained in precision measurements at LEP and in other experiments ||. This prompts us to 
test the approach to the continuum limit of QCD thermodynamics by testing the constancy of T c /A-^, where A-jjg- 
is the QCD scale parameter extracted in the MS scheme. 

In the limit when all quark masses are zero (or infinite), QCD has only one dimcnsionlcss parameter the coupling, 
a s — g 2 /47T. Quantum corrections transmute it into a momentum scale. This scale is either given explicitly as the 
QCD parameter A, or specified implicitly, as the value of the coupling a s (/i) at scale /j,. 

In the lattice regularisation of QCD, the value of the lattice spacing (a) is determined by the bare coupling (3 = 6/g 2 - 
However, 6/(3 is not a good expansion parameter for perturbation theory. It is more appropriate to define the coupling 
through some physically motivated process. In one definition, called the V-scheme H, a s at scale 3.4018/et is found 
from the logarithm of the plaquette value (V) Q through the formula 

47T 

- ]nV = —a v (1 - (1.1897 + 0.071N f )a v ) . (1) 
In another definition, the E-scheme Q , the renormalized coupling at scale a is 

"-5(1-"). < 2 » 

We choose to work in these two schemes, as well as the MS scheme whose relation to the V-scheme has been worked 
out H . At 2-loop order we can write 

aA = kR (l/47r/? av) , where R 2 (x) = exp(-x)x l3l/f3 ° . (3) 

The constant k depends on the scheme, being unity in the E-scheme and 3.4018 in the V-scheme. The function R is 
obtained by integrating the two-loop beta function, (3(g) = —(3og 3 — (3ig 5 . 

The QCD phase transition temperature is determined by tuning the bare coupling on lattices with N t <C N s (where 
N t is the number of sites in the Euclidean time direction and N s that in the spatial directions). Then, T c — l/a c N t , 
where a c is the lattice spacing at the coupling where the phase transition occurs. Our strategy is to compute the 
renormalised coupling a v at these bare couplings and hence determine Tj A— using eq. ([|) , and the known ratio of 
A-jjg and the scale A in any other scheme. 

In quenched QCD with the usual Wilson action, the critical bare couplings, (3 C , have been determined for 2 < iV f < 16 
H . The main source of systematic uncertainty in the older data arises from the fact that the thermodynamic limit 
N s — ► oo was not taken. Later data have taken this limit, and we only used these to study scaling. The statistical 
errors in these later studies are also much smaller, and hence they are able to test the scaling hypothesis much more 
stringently. 
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1 V is defined to be one third of the real part of the trace of the product of four link matrices taken in order around a plaquette. 
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We extracted a s from the plaquette values listed in using eqns. (|l],|j). Values of hiV at C were obtained by cubic 
spline interpolation. Statistical errors in the interpolated values were found by propagation. We probed systematic 
errors in the interpolation by the change in a s on removal of some of the knot points. The results of our analysis are 
shown in Figure [l]. T c /A— is constant down to (3 C for JV* = 3. A large part of the error in this ratio comes from the 
errors in the measurement of V, which, while small, are exponentiated in A. 

More detailed results are shown in the right hand panel of Figure [|. Small scaling violations at these couplings, 
seen in measurements at T — 0, have been attributed to finite lattice spacing errors JTo[ ] . On replacing the scaling 
function R in eq. (||) by R(a v ) = R{a v ) [l 4- C2<2 2 + 04a 4 ] , where a = R(a v )/R(a v ) and a v is determined at (3 = 6, 
it was found jll[] that data on a^fa could be quantitatively described down to = 5.4. As shown in the figure, in 
the MS-scheme we can fit a constant value of T c /Ajjg using the critical couplings determined for 3 < N t < 14 with 
good confidence — \ 2 = 0.6 for 4 degrees of freedom. The statistical error in this estimate is about 1%. Taking into 
account the variation between values obtained in different schemes, and for fits with N t > 3 or 4, we quote 

T c = (1.15 ± 0.05)A— (Quenched QCD). (4) 

This error estimate now includes not only the statistical errors but also the systematic uncertainties above. Since 
experiments cannot access a quenched theory, the value of A-j^ in MeV units is hard to pin down. When comparing 
different estimates, in MeV units, of the same scale in quenched QCD this universal problem should be borne in mind. 
Here we convert T c into physical units using the value A^g « 250 MeV as a mean of the values obtained in different 
schemes in |nj by assuming yfa = 465 MeV. This gives T c ss 285 ± 10 GeV, in reasonable agreement with other 
estimates, such as that found in a recent study using an RG-improved action |12j. 
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FIG. 1. The panel on the left shows T c /Aj^ as a function of N t in the V-scheme. For Nf = we show tests of scaling for 
both the bare lattice coupling and a v . For the former the errors come only from the determination of the critical coupling, 
whereas for the latter they also include errors in the determination of V . For Nf = 2 with m/T c — 0.1 and Nf = 4 with 
m/T c — 0.08 only the latter is shown. The panel on the right shows details of T c /Ajjs for Nf — in the MS-scheme. 



This analysis leads us to investigate the extraction of thermodynamic quantities in the QCD plasma using renor- 
malised perturbation theory. Failure of scaling would then be a direct signal for lattice artifacts such as power 
corrections in a. We examine the energy density (E) and the pressure (P). These can be written in terms of the 
difference Aj = Vi — Vq between the spatial (i — s) and temporal (i = t) plaquettes and their zero temperature 
counterpart Vq. E is defined by the formula 



E A 

2* = ™ C N? 



47r«c 



«A S + c' t A t ) 



(5) 
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The anisotropy coefficients c' s and c' t are known to 1-loop order g jl3| . A measure of deviations from ideal gas behaviour 
is 

F — 3P 

A = = l2N c Nt{c' s + c' t )(A t + A.), (6) 

which can be combined with eq. (jsj) to give the pressure. The sum rule <? 3 (cl, + c' t ) = /3(g) jl3|], allows us to evaluate 
A beyond 1-loop order. In fact, part of the finite a corrections can be incorporated into A by evaluating (3(g) using 
R instead of R. However, we avoid this approach, since we need to use c' s and c' t consistently in eqs. (pIS. In the 
following we use only the 1-loop results for these coefficients and work in the MS scheme. 




FIG. 2. (E~3P)/T 4 for a QCD plasma evaluated as a function of T/T c at different lattice spacings when a-jj^ sets the scale. 
The continuum curve is not continued into the region where finite volume effects may be important. Errors on the continuum 
extrapolation are roughly similiar to those at N t = 8. 



Raw data on Vo, s ,t at a range of couplings measured on Nt X (4/Vt) 3 lattices (Nt — 4, 6 and 8) are tabulated in 
||. We have used these to evaluate thermodynamic quantities only for T < 4T C , because the finite spatial volumes 
of the lattices used may cause spatial deconfinement above 4T C . As expected, A varies as 1/N 2 at fixed T, i.e., as 
a 2 . This power correction can be removed through the use of improved actions fbl . A shows a peak at T w l.lTc 
in the continuum limit. However, the location of the peak is uncertain because there is no unique definition of T c 
at finite volumes. Different definitions, which all coincide in the thermodynamic limit, give different values of the 
pseudo-critical point on finite volume systems p5[ . Using coarse lattices, we have estimated that in SU(3) theory for 
N s = ANt this inherent uncertainty in the critical coupling may be as much as 8(3 ~ 0.005-0.01: much larger than the 
statistical error for any given definition of the critical point. For similar reasons the value of A at the peak cannot 
be reliably extracted without taking into account finite volume effects. In the range 1.3T C < T < 4T C , finite volume 
corrections are expected to be small. For such T, A is monotonically decreasing. However, many quantities of interest 
can be extracted only near T c . The results of a finite size scaling study which does this will be reported elsewhere. 

We have also evaluated E/T 4 and found that it scales to the continuum limit as a 2 — l/N 2 . In the temperature 
range where finite volume effects are strong, the large value of A gives a negative value to P evaluated using the 
formulae in eqs. (^,^|). As we have argued already, an evaluation of thermodynamics in this range of T requires better 
control over finite volume effects. However, the problem of negative pressures is also avoided if the continuum limit 
of P is found by combining E/T 4 and A in the limit a = 0. This gives nearly vanishing P in the region of the peak 
in A. 



2 At this order the coefficients remain unchanged in going from the lattice scheme, g 2 — G/(3, in which they have been 
computed, to the MS or V-scheme, but change in the E-scheme. In the absence of a two-loop computation, we have evaluated 
the renormalised couplings at the scale appropriate to the plaquette. 
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FIG. 3. E/T 4 and 3P/T 4 for a QCD plasma evaluated as explained in the text, using the data of ||. The dashed curves are 
smooth cubic spline fits. The dotted horizontal line is the ideal gas result in the continuum limit. 

In Figure || we show our estimates of the continuum limit of E and P in the range 1.3T C < T < 4T C , where finite 
volume effects are small. At 4T C our estimate of the energy density is about 15% lower (and the pressure about 19% 
lower) than for an ideal gas. At this temperature, our estimates of E and P are about 5% lower than the earlier 
determination which used the integral method ||. While statistically insignificant at present, this difference is due 
entirely to the treatment of the region around the peak of A [J The speed of sound in the QCD plasma (c s ) is slightly 
lower than the ideal gas value (c 2 s — 1/3) at 4T C (but consistent with it within errors), and falls to c 2 s w 0.1 ± 0.025, 
at 1.5T C . 

Our attempt to measure the pressure using perturbativcly determined coefficients may seem mysterious when 
current practice has converged on the use of an "integral" method ||. As we have mentioned, these two methods 
agree within error bars for T > 2T C but differ near T c . The integral method, by construction, gives a continuous 
pressure across T c . However, if true, this is important information and needs independent verification. The latent 
heat density at the phase transition, AQ/T^ = AS/T^, where AS is the jump in the entropy density at T c . In general 
this is not the same as the change in the energy or enthalpy density at the phase transition. These contain an extra 
term due to the change in specific volume at fixed pressure or the change in pressure at fixed volume. In a quenched 
lattice simulation, the change in volume is negligible when changing T by an infinitesimal amount across T c . As a 
result, the pressure should change discontinuously at T c . This is analogous to boiling water in a closed vessel — the 
pressure jumps at the phase transition. Such a jump in pressure is related to the change in specific volume, and 
hence the slope of the phase boundary in the (P, T) plane through the Clausius-Clapeyron equation. This physics is 
missed when the integral method is used to extract the pressure. This issue is currently being investigated in a lattice 
simulation whose results will be reported elsewhere. 

While the observation of scaling on coarse lattices is interesting in quenched QCD, the real pay off would be in the 
study of finite temperature QCD with dynamical light quarks. Since there are inherent difficulties with simulating 
QCD with dynamical massless quarks, all lattice studies have used quarks with mass m > 0. In such theories, scaling 
should be tested at fixed m/Aj^ (or equivalently, fixed ratio of m and any hadronic mass scale). 

Measurements of (3 C in QCD with 4 flavours of dynamical staggered quarks have been performed for 2Vj = 4, 6 |l(t | 
and 8 p7[ . This last measurement was done with m/T c = 0.08. The simulations at smaller N t were done at several 
values of m, enabling us to find (3 C at m/T c = 0.08 by interpolation. Plaquette values were taken from a recent finite 
temperature simulation at fixed m/T c Gq|. Since it is known that A t , s /Po are less than 0.1%, even near T c , we have 
used (P s + Pt)/2 to determine a s . The difference of cv s measured using V s and Pt is taken as an estimate of its error. 
We used eq. (||) to set the scale. T c /Aj^, when determined through the bare coupling, changes from 8.2 to 4.7 in 
going from Nt = 4 to 8. As shown in Figure El when a v is used, 2-loop scaling works much better. 




3 However, note that out use of the 1-loop expressions for c' 3 and c' t , and the consequent inability to make an optimal choice 
of scale is the largest uncertainty in this method for T > 1.5T C . Such uncertainties can be reduced by appropriate 2-loop 
computations. 
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The phase transition in 2-flavour QCD has been studied in greater detail (see M for a recent compilation of data). 
Measurements of (3 C using dynamical staggered quarks exist for N t < 12 Earlier simulations with Wilson 

quarks (which have order a lattice artifacts) showed that pion masses were rather high compared to those obtained 
with staggered quarks. This problem becomes less acute on using improved actions for Wilson quarks, and finite 
temperature simulations have now been performed with such improved actions pT| , p2"| . There have also been some 
studies with domain wall Fcrmions. 

We have fixed a v and set the scale using published plaquette values for staggered quarks at several bare quark 
masses pp| , p3| . T c /A— for Nf = 2 shown in Figure [l] are based on the subset of the data which uses staggered quarks 
at m/T c = 0.1. In contrast to the near constancy of this ratio as shown in the figure, T c /Ajf^ computed from the 
same data using the bare lattice definition of g falls from 3.9 to 2.8 in going from N t — 4 to 8. To extend this test to 
other values of m/T c we have to interpolate between plaquette values for various quark masses. An upper limit for 
the error in this procedure is the actual change in V between extreme values of the quark masses at which they are 
measured. This varies between a few parts in a thousand at (3 = 5.26 to about 2% at f3 — 6. This uncertainty in V 
translates into a similar magnitude of uncertainty in T c j 'A-jjg and is much smaller than the change when using 6//3 
across a similar range of quark masses. Excellent scaling of T c /Ajys 1S seen also for other values of m/T c . 

In order to obtain a physically relevant value of T c , it is necessary to extrapolate T c / Ajrs to measured values of the 
hadron masses. It would be most interesting to perform this extrapolation in the quark or pion masses. However, 
this needs control over the critical exponents of the theory — a task we do not attempt here. Instead we choose to 
extrapolate T c to the physical region in terms of M p There are two reasons for this. First, the ratio T c /M p is 
known to be nearly constant. Secondly,M p is quite sensitive to finite lattice spacing effects. The linearity of the plot 
of T c / ' Ajjg against M p / k-jfg in Figure [| then indicates that, for the chosen data set, finite lattice spacing effects in 
M p are under reasonable control. The value of A— also needs to be specified. This depends on how many active 
flavours are present at the scale under consideration. The world average of Ai^L (at scales high enough for 5 active 

flavours) is 219+23 @- At lower scales, with only three active flavours, Ai^l = 343+28: using the prescription of jj| to 
match across flavour thresholds. 
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FIG. 4. T c /h-jjs plotted against M p /A-^ and extrapolated to the physical value. The band shows the 1-er errors on the 
best fit line to the measurements. 

The physical value of T c is obtained by extrapolating to the real- world value of M p /A^L. This is done by fitting a 
straight line to the data including staggered and improved Wilson quarks Q . We have taken the two sets of data 
together since they both have cutoff effects of order a 2 . The extrapolation gives 

T c /A— = 0.487 ± 0.023, and T c = 167 ± 9±\l (N f = 2 QCD). (7) 

The error in the ratio above reflects only the statistical errors in various measurements. The first error in T c is purely 
from extrapolation and the second set from the errors on Since they come from independent sources, it is 

possible to add them in quadrature. This result is consistent with a recent estimate using 0(4) critical indices to scale 
T c as the pion mass is taken to its physical value [^2|. However, there are possibly large systematic uncertainties. 
Currently the least well-understood problem is whether extrapolating one hadron mass to its measured value also 
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takes all other hadron masses to their correct values. Until this issue is settled, all estimates of T c must be considered 
preliminary. 




FIG. 5. a s at 2-loop order in the Nf = and Nf = 2 theories at the scale optimised for the three-dimensional gauge 
coupling in dimensionally reduced theories, shown as a function of the temperature. The bands arise from the quoted errors 
in the measurement of T c /Ajj-g. This plot is used to limit the validity of dimensional reduction for Nf = 2 as explained in the 
text. 



We end with a small application of our measurements of T c / h-^. This addresses the question of the temperature 
range where dimensional reduction (DR) is expected to be valid. Perturbative matching of the T > four dimensional 



theory with DR fixes the effective couplings in the latter 
gauge theory agrees with its DR version for T > 2T C [p5 



!4|. There is good numerical evidence that the SU(3) pure 
If the value of the strong coupling were the determining 



factor in the agreement, then this would imply that for Nf = 0, a s is small enough at 2T C to trust this expansion. If 
for Nf = 2 no other physics becomes more important, then the matching should be equally reliable when the coupling 
is as small. Now we ask for the temperature at which a s in the Nf = 2 theory takes on the value that it had at 2T C for 
Nf = 0. Using the values of T c /Ajy^ already extracted, we are in a position to do this. The optimal scale choice for 
a s has been investigated in p4[ . This depends, of course, on Nf and the quantity for which the perturbation theory 
is optimised. Choosing the scale which optimises the three-dimensional gauge coupling we find the dependence of a s 
on T shown in Figure |5[ This tells us that for Nf = 2 DR can be used only for T > 6.5T C . Using scales optimised 
for other quantities, and varying the computation of the strong coupling from 2-loop to 3-loop order, we find that the 
lower limit of the range of validity of DR for Nf — 2 varies between 6T C and 8T C . 

In summary, we demonstrated that the lattice data on strong interaction thermodynamics obey QCD scaling 
relations very well and allow some continuum physics to be extracted on fairly coarse lattices. T c can be measured 
with precision of about 4%, taking into account both statistical errors as well as systematic uncertaintites such as 
renormalisation scheme dependence. Simulations including dynamical quarks, when extrapolated to physical values 
of the p meson mass give statistical errors of about 4% in T c j ' A-^-g. Renormalised perturbation theory at 1-loop order 
seems to be able to give a good account of the pressure, energy density and the speed of sound for T > 1.5T C , leading 
to hopes that improved actions and higher loop orders can yield continuum physics at small expense. A study of finite 
volume effects on thermodynamics closer to T c will be reported elsewhere. 

It is a pleasure to thank Urs Heller for making available the full results of the fits reported in jll], and Rajiv Gavai 
for many discussions. 
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